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Abstract. 

We discuss the gravitational self-force on a particle in a black hole space-time. For 
a point particle, the full (bare) self-force diverges. The metric perturbation induced by 
a particle can be divided into two parts, the direct part (or the S part) and the tail part 
(or the R part), in the harmonic gauge, and the regularized self- force is derived from 
the R part which is regular and satisfies the source-free perturbed Einstein equations. 
But this formulation is abstract, so when we apply to black hole-particle systems, 
there are many problems to be overcome in order to derive a concrete self-force. These 
problems are roughly divided into two parts. They are the problem of regularizing the 
divergent self- force, i.e., "subtraction problem" and the problem of the singularity in 
gauge transformation, i.e., "gauge problem". In this paper, we discuss these problems 
in the Schwarzschild background and report some recent progress. 
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1. Introduction 

Thanks to recent advances in technology, an era of gravitational wave astronomy has 
almost arrived. There are already several large-scale laser interferometric gravitational 
wave detectors that are in operation in the world. Among them are TAMA300 p], 
LIGO [2j, GEO-600 [SI, and VIRGO 4j is expected to start its operation soon. The 
primary targets of these ground-based detectors are inspiralling compact binaries, which 
are expected to be detected in the near future. 

There are also space-based interferometric detector projects. LISA is now on its 
R&D stage 0, and there is a future plan called DECIGO/BBO ^[Zj. These space- 
based detectors can detect gravitational waves from solar-mass compact objects orbiting 
supermassive black holes. To extract out physical information of such binary systems, 
it is essential to know the theoretical gravitational waveforms with sufficient accuracy. 
The black hole perturbation approach is most suited for this purpose. In this approach, 
one considers gravitational waves emitted by a point particle that represents a compact 
object orbiting a black hole, assuming the mass of the particle {^) is much less than 
that of the black hole (M); fi<^M. 

In the lowest order in the mass ratio (fi/M)^, the orbit of the particle is a geodesic 
on the background geometry of a black hole. Already in this lowest order, by combining 
with the assumption of adiabatic orbital evolution, this approach has been proved 
to be very powerful for evaluating general relativistic corrections to the gravitational 
waveforms, even for neutron star-neutron star (NS-NS) binaries where the assumption 
of this approach is maximally violated [S]. 

In the next order, the orbit deviates from the geodesic because the spacetime is 
perturbed by the particle. We can interpret this deviation as the effect of the self-force 
of the particle on itself. Since it is essential to take account of this deviation to predict 
the orbital evolution accurately, we have to derive the equation of motion that includes 
the self-force. The gravitational self-force is, however, not easily obtainable. There are 
two main obstacles. 

First, the full (bare) metric perturbation diverges at the location of the particle, 
which is point-like, hence so does the self-force. The full metric perturbation can 
be formally divided into two parts; the direct part that comes from the perturbation 
propagating along the light cone of the background spacetime directly from the source 
particle, and the tail part that arises from the curvature scattering of the perturbation, 
which exists within the light cone. The self-force is given by the tail part of the metric 
perturbation which is regular at the location of the particle jHl QHI- Recently, this is 
reformulated in a more elaborated way by Detweiler and Whiting ^T]. In this new 
formulation, the direct part is replaced by the S part and the tail part is replaced 
by the R part. The R part is equivalent to the tail part as far as the self- force 
calculation is concerned, but which has a much better property that it is a solution 
of source-free linearized Einstein equations. Thus, one has to identify the R part of the 
metric perturbation to obtain a meaningful self-force. However, by construction, the 
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R part cannot be determined locally but depends on the whole history of the particle. 
Therefore, one usually identifies the direct part (or the S part) which can be evaluated 
locally in the vicinity of the particle to a necessary order, i.e., by local analysis, and 
subtract it from the full metric perturbation. This identification of the S part is the 
subtraction problem. 

Second, the regularized self-force is formally defined only in the harmonic (Lorentz) 
gauge, in which linearized Einstein equations are hyperbolic, which enables us to define 
the S part and the R part uniquely. On the other hand, the metric perturbation of a black 
hole geometry can be calculated only in the radiation gauge in the Kerr background |12j . 
or in the Regge- Wheeler (-Zerilli) gauge in the Schwarzschild background fHlE]- Hence, 
one has to find a gauge transformation to express the full metric perturbation and the 
S part in a same gauge, to identify the R part in that gauge correctly. This is the gauge 
problem. 

As a method to solve the subtraction problem, extending earlier work ^T^l ITH] , 
Mino, Nakano and Sasaki formulated two types of regularization methods, namely, 
the 'power expansion regularization' and the 'mode- decomposition regularization' ^Tj 
HHj. In particular, the mode-decomposition regularization was found to be quite 
powerful. Meanwhile, Barack and Ori independently developed the 'mode-sum 
regularization' ^3 1201 1^- Although the mode-decomposition regularization and the 
mode-sum regularization were quite different in their formulations, they are essentially 
the same in the sense that the spherical harmonic expansion is used to regularize the 
divergent metric. In fact, it was shown that the two methods give the same result for 
the S part in the harmonic gauge, in the form of a (divergent) spherical harmonic series, 
for an arbitrary orbit in the Schwarzschild background |22j . 

However, it seems extremely difficult to solve the metric perturbation under the 
harmonic gauge because the metric components couple to each other in a complicated 
way [23 . This is one of the reasons why the gauge problem is difficult to solve. Recently, 
Barack and Ori [21] gave a useful insight into the gauge problem. They proposed an 
intermediate gauge approach in which only the direct part of the metric in the harmonic 
gauge is subtracted from the full metric perturbation in the Regge-Wheeler (RW) gauge. 
They then argued that the gauge- dependence of the self-force is unimportant when 
averaged over a sufficiently long lapse of time. Using this approach, the gravitational 
self-force for an orbit plunging into a Schwarzschild black hole was calculated by Barack 
and Lousto [25]. 

In the case of the Kerr background, which is of our ultimate interest, there is 
additionally a much more intricate technical issue of how to treat the spheroidal 
harmonics: Neither the eigenfunction nor the eigenvalue has simple analytical 
expressions, and they are entangled with the frequency eigenvalues of Fourier modes. 
Further, the only known gauge in which the metric perturbation can be evaluated 
is the radiation gauge formulated by Chrzanowski J2]- However, the Chrzanowski 
construction of the metric perturbation becomes ill-defined in the neighborhood of the 
particle, i.e., the Einstein equations are not satisfied there [21]. Some progress was 
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made by Ori [^Hl to obtain the correct, full metric perturbation in the Kerr background. 
The regularization parameters in the mode-sum regularization for the Kerr case are 
calculated by Barack and Ori [77]. Although our ultimate goal is to overcome these 
difficulties altogether, since each one is sufficiently involved, we choose to proceed step 
by step, and focus on the Schwarzschild background. 

In this paper, in order to obtain the gravitational self-force, we consider a particle 
orbiting a Schwarzschild black hole, and propose a method to calculate the regularized 
self-force by solving the subtraction and gauge problems simultaneously. Namely, 
we develop a method to regularize the self-force in the Regge- Wheeler gauge. The 
regularization is done by the mode-decomposition regularization, and we introduce the 
concept of a "finite gauge transformation" of the S part j2Hl 1211 • It is noted that, 
although our approach is philosophically quite different from the intermediate gauge 
approach by Barack and Or \^A^ , practically both approaches turn out to give the same 
result as far as the calculation of the gauge transformation is concerned. 

The paper is organized as follows. In Sectional we review the mode-decomposition 
regularization. We first derive the direct part of the self-force (which is equivalent to the 
S part in the coincidence limit of a field point with a point on the particle trajectory) 
by local analysis. Then we transform it to an infinite series of spherical harmonic modes 
and derive the regularization counter terms for a general orbit. In Section El further 
focusing on a circular orbit, we perform a gauge transformation from the harmonic 
gauge to the RW gauge, and present the first post-Newtonian order self-force. One of 
the reasons why we had to focus on circular orbits is that it is formidable to obtain a 
closed analytic expression of the full metric perturbation for a general orbit even under 
the post-Newtonian expansion. In Section IH turning back to a simpler case of the 
self-force of a scalar-charged particle, we present a new regularization method that can 
apply to an arbitrary orbit. This new method is a variation (or an improved version) 
of the mode-decomposition regularization. The essential point is to divide the full field 
in the frequency domain into two parts, the S part and R part, in such a way that 
the divergent part is totally contained in the S part whose frequency integral can be 
analytically performed. 

2. Mode-decomposition regularization 

In this section, we review the mode-decomposition regularization method, and derive 
the regularization counter terms for an general orbit in the Schwarzschild background. 
For simplicity, we consider a scalar-charged particle which sources a scalar field ip. The 
extension to the gravitational case is straightforward. The only complication is the 
presence of additional tensor indices associated with the metric perturbation. 

Since the tail part depends non-locally on the geometry of a background spacetime, 
it is almost impossible to calculate it directly. However, for a certain class of spacetimes 
such as Schwarzschild/Kerr geometries, there is a way to calculate the full field generated 
by a point source. Considering a field point shghtly off the particle trajectory, it is then 
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possible to obtain the tail part by subtracting the locally given divergent direct part 
from the full field. 

Let X = z{r) be a trajectory of the particle with r beging the proper time along 
the trajectory Then the field (p{x) generated by the particle is given by 

-q j dTG^\x,z{T)), (1) 

where G^^^^{x,x') is the retarded Green function satisfying the Klein-Gordon equation, 

V^V.G^^^x, x') = - ^^'^^^''') . (2) 

The self-force at a; = z(ro) is schematically given by 

F,(to)= hin (3) 

I— >z(to) 

where 

F^[^'^\x) = F^[^'^\x) - , {x ^ z{t)) . 

The symbols v^^"", v^'^''^ and stand for the full field, the direct part and the tail part, 
respectively. | Both ip^"^^ and (z?*^"' diverge in the coincidence limit x — > z{t), while ip^^""^ 
is finite in the coincidence limit. is a tensor operator on the field, and is defined 

as 

FM=<lPa^^P^, (4) 

where Pi/ = S^ + VaV^ is the projection tensor with V"' being an appropriate extension 
of the four velocity v'^{tq) off the orbital point. 

In what follows, we use the Boyer-Lindquist coordinates for the Schwarzschild 
background, 

ds^ = - (l - V") "^^^ + ~ V") ' + ""^ ("^^^ + ''''' ' 

and denote a field point by x = {t, r, 9, 4>} and a point on the orbit by z{tq) = zq — 
{to, To, 6*0, (po}. 



2.1. Basis 

Even in the Newtonian limit, the integration of an orbit involves an elliptic function. 
Thus, numerical computations will be necessary at some stage of deriving the self-force. 
However, it will not be easy to perform the regularization of the divergence numerically. 
The idea to overcome this difficulty is to replace the divergence by an infinite series, 
each term of which is finite and analytically calculable. 

The derivation of the direct part (p^^'^{x) is one of the main issues in the 
regularization calculation. The direct part of the scalar field is obtained by integrating 

I In this section, wc consider the direct part instead of the S part, because the direct part has a sUghtly 
simpler form for local analysis. 
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the direct part of the retarded Green function with the scalar charge. The direct part 
of the retarded Green function C^^^ is given in a covariant manner as jHO] 

G^"(x,x')= - ^d[i:{x),x'y/\{x,x')5{a{x,x')), (6) 
47r 

where a{x,x') is the bi-scalar of half the squared geodesic distance, A(x, x') is the 
generalized van Vleck-Morette determinant, S(a;) is an arbitrary spacelike hypersurface 
containing x, and 6'[S(x),x'] = 1 — 6[x' , is equal to 1 when x' lies in the past of 

and vanishes when x' lies in the future. The physical meaning of the direct part 
is understood by the factor 0\T,{x),x']5{a{x,x')) in Eq. (jH)). Since cr(x, x') describes the 
geodesic distance between x and x', the direct part of the Green function becomes non- 
zero only when x' lies on the past lightcone of x. Hence the direct part describes the 
effect of the waves propagated directly from x' to x without scattered by the background 
curvature. 

For the actual evaluation of the direct part, a couple of methods have been proposed. 
In Refs. jHH the direct part of the field is calculated by picking up a limiting 
contribution in the full Green function from the light cone as 

/■oo 

¥P^"(x)=lim/ dTG^"^\x,z{T))S{T), (7) 

where G^"^^ is the retarded Green function, S{t) is the scalar charge density, and Tret (a;) 
is the retarded time defined by the past light cone condition of the field point x as 

e[S(x), z(rret)]5(a(x, z{r,,,))) = . (8) 

However, the calculation seems rather cumbersome when we apply this method to a 
general orbit. 

In Ref. pn], the direct part was evaluated using the local bi-tensor expansion 
technique. Using the bi-tensor, the direct part is expanded around the particle location 

as 

+ 0{y') , (9) 

where the letters /i, ■ ■ ■ are used for the indices of the field point x, a, /5, ■ ■ ■ for 
the indices of the orbital point and f"(r) is the orbital four velocity at z{t). The 
order of the local expansion is represented by the powers of ?/, where y is the coordinate 
difference between the field point x and the orbital point Zq] = x^ — Zq. Because the 
full force is quadratically divergent, we must carry out the local bi-tensor expansion of 
the full field up through Oiy). By evaluating the local coordinate values of the relevant 
bi-tensors, we obtain the local expansion of the full force in a given coordinate system. 
As described in Ref. ^H], this may be done in a systematic manner, and it is possible 
to obtain the explicit form of the divergence for a general orbit. However, the problem 
is how to decompose it into an appropriate infinite series. 



(X,^(ri.et))w°(rret) 
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2.2. Mode- decomposition regularization 

The mode-decomposition regularization method is based on the spherical harmonic 
series expansion ^H]- There is a delicate problem in this approach. The exact 
decomposition calculation usually needs the global analytic structure of the field so that 
we can uniquely define each term in the infinite series. On the other hand, the divergent 
direct part to be subtracted is given only by local analysis of the field El EOl 
Thus the direct part is defined (and known) only in a local neighborhood of the particle. 
Because of this, to obtain a harmonic series expression of the direct part, we need a global 
extension of the direct part. But this extension cannot be done uniquely in practice. 
Nevertheless, the final result of the self-force should be unique. Although we have no 
explicit proof for the uniqueness of the regularization counter terms for the self-force, 
the uniqueness is intuitively apparent because the extension is performed in such a 
way that the local behavior is retained with sufficiently accuracy. In addition, it is 
reassuring to find that the resulting counter terms, presented at the end of this section, 
agree completely with those obtained in the mode-sum regularization by Barack and 

ori iinii2niinii22i- 

The harmonic decomposition is defined on a two-sphere. However both the direct 
field and the full field have divergence on the sphere containing the particle location, the 
mode decomposition is ill-defined on that sphere. Therefore, we perform the harmonic 
decomposition of the direct and full fields on a sphere which does not include, but 
sufficiently close to the orbit. The steps in the mode-decomposition regularization are 
as follows. 

1) We evaluate both the full field and the direct field at 

x = {t,r,^,0}, (10) 
where we do not take the coincidence limit of either t or r 

2) We decompose the full force and direct force into infinite harmonic series as 

Im 

where F„ [(^f""/'^''"] (x) are expanded in terms of the spherical harmonics Y(m{6,4>) 
with the coefficients dependent on t and r. For the direct part, the harmonic 
expansion is done by extending the locally defined direct force over to the whole 
two-sphere in a way that correctly reproduces the divergent behavior around the 
orbital point zq up to the finite term. 

3) We subtract the direct part from the full part in each i, m mode to obtain 

^^b""] = (^^b'-^i - • (13) 

Then, we take the coincidence limit x zq. Here we note that one can exchange 
the order of the procedure, i.e., first take the coincidence limit and then subtract, 
provided the mode coefficients are finite in the coincidence limit. 
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4) Finally, by taking the sum over the modes, we obtain the self-force, 

K{r,)=Y^F':-[^'^\z,). (14) 



2.3. Decomposition of the direct part 

The advantage of using the expression Q for the direct part is that we have a systematic 
method for evaluating it. Here we describe a method to evaluate the bi-tensors in a 
general regular coordinate system. 

Before we consider the local expansion in a given coordinate system, we calculate 
the derivative of Eq. Q, and derive the direct part of the force with the local bi-tensor 
expansion using the equal-time condition. 







^a(x,.(r)) 



(15) 



T=Taq(x) 

We define an extension of the four-velocity off the orbit by 

dz' 



V°(X) := ^as(a:,2;eqXq; ^eq = 2:(req(x)) , v"^^ = — — , (16) 



where Qaa is the parallel displacement bi-vector. 
Using the formulas in Refs. [HI ED] , we have 



+ \^'^Ra^^-5{Zeq)vlqCr''^{x, 2;eq)t^fq \ + 0{y) , (17) 



2a(a;,2eq), (18) 



K 



— Cr^^(x, Zeq)^"q^eq 



= 1 + -R^-p^sizeciK^a'f^ix, Zeq)<a'^(x, z,^) + 0{y^) . (19) 

The bi-tensors necessary for the evaluation of the direct force ()17|) are (t(x, x) and 
gaa{x,x), for which we consider the local coordinate expansion of these bi-tensors 
around the coincidence limit x ^ x. In a general regular coordinate system, a{x, x) and 
gaa{x,x) can be expanded as 



^^?a/3(x)l/"^+ V ^A,i,2...,„(x)y-'-'-"", (20) 
n=3,4,--- 

+ -,B^^\P-P-...p.{x)y^'^'---^" ■ (21) 



n=3.4 



n=l,2,--- 



° - = (x" )(x° -x" )■■■. (22) 

To calculate the reaction force to a particle, it is enough to know the expansion 
coefficients of n = 3,4 of Eq. ()20p and n = 1,2 of Eq. (PT|) . These are shown in 
Appendix A of Ref . [IH] . The local expansion of the force (|T7jl on the Boyer-Lindquist 
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coordinates is quite tedious, though systematic. However, most of the terms give no 
contribution to the harmonic coefficients in the coincidence hmit x — > Zq- Below, we 
shall focus on the terms that are non- vanishing in the coincidence limit. 

Without loss of generality, we may assume that the particle is located at = 7r/2, 
00 = at time to- Since the full force is calculated in the form of the Fourier-harmonic 
expansion and the Fourier modes are independent of the spherical harmonics, we may 
take the field point to lie on the hypersurface t = to in the full force. Hence we may 
take t = to before we perform the local coordinate expansion of the direct force. That 
is, we consider the local coordinate expansion of the direct force at a point {to,r, 6',0} 
near the particle localtion {to, tq, Tr/2, 0}. 

The local expansion of the direct force in the Boyer-Lindquist coordinates can be 
done in such a way that it consists of terms of the form. 



^2n4+l ' 

b 



(23) 



^:=V2rJa-cose + ^{(j)-(j)'y] , (24) 

R:=r-ro, Q := 9 - ^ , (25) 

where ni, n2, n^, are non- negative integers, and a, h and defined by 

a := 1 + /° o 7 ^° (26) 

^ - ' (27) 



„2~r~;^"^-^' (2^ 



where £ := —gudt/dr, C := g^jj^dcp/dr, and Ur := grrdr/dr, and 9 is the relative angle 
between (9,4)) and (7r/2,0'). 

There are two apparently different terms in the covariant form of the direct force 
given by Eq. (jl7|) : the ffist term in the curly brackets exhibiting the quadratic divergence, 
and the second term proportional to the curvature tensor that appears to be finite in the 
coincidence limit. In the local coordinate expansion, the second term will give terms of 
the form R/^ or 0/^. As discussed in Section the harmonic coefficients of R/^ vanish 
in the coincidence limit, while those of (p/^ are finite but they give no contribution to 
the final result when the infinite harmonic modes are summed up after the coincidence 
limit is taken. Hence we may focus on the ffist term in the curly brackets of Eq. (fTTj) . 

Since the orbit always remains on the equatorial plane, the force is symmetric under 
the transformation 9 ^ 7i — 9, which implies there is no term proportional to odd powers 
of G. Hence we only need to consider the case of n2 being an even number in the general 
form given by Eq. Then the factor B"^ may be eliminated by expressing in 

terms of ^, R and 0, and we are left with terms of the form, 

(29) 



2n4+l 
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Explicitly, we find 



1 (ro - 2M)Sur 1 ^ 2(ro - 2M)SC'^Ur 



2 ri e r^o e 

3(ro -2M)££V </^^~ 



2 e 
dir^ / R rjS' R I 

^ ^Uo(ro-2M)e (ro-2M)2e ""^^ 



2 rg ^ 2ro-2M^ 



+ 



2 e 2 ^ 



(30) 



2 rj] ^3 ro-2M^3 

3 (rg + £^)£^ 0^ 3 g^/:^ 0^ 

2 2ro-2Me' 
= , (32) 

Ff = g(^-£«.|-(ro^ + £^)^ 

1 (ro - 2M)/:m, 1 1 (ro - 2M) (r^ + 4£2)£«^ 02 



, 3 (ro-2M)(rg + £^)£V 0^^ 

+ 2 ro^ ej' ^ ^ 

where F^" = . The absence of F^" is because of the symmetry of the background. 



2.4- Regularization counter terms 

What we have to do now is to perform the harmonic decomposition of the components of 
the direct force given above. To do so, we note the following important fact. Apart from 
the trivial multiplicative factor of R^^ which is independent of the spherical coordinates, 
the terms to be expanded in the spherical harmonics are of the form 0"3y^2n4+i^ qj. 
{(f) — 0')"'3/^^"*+^. To the order of accuracy we need, the factor {(f) — (f)'Y^ may be 
eliminated by replacing it to an equivalent 0-derivativc operator of degree acting on 
^2n3-2n4-i^ which is further converted to a polynomial in m after the harmonic expansion 
Qf ^2n3-2n4-i_ thc Only basic formula we need is the harmonic expansion of ^^p-i 

where p is an integer. Note that, apart from the term 6(0 — 0')^/2 in ^ with respect to 
which we expand ^ in a convergent infinite series, ^^^"^ is defined over the whole sphere 
to allow the straightforward harmonic decomposition. The result to the leading order 
in the coincidence limit a — > 1 + is 

[-^J - (a-co.. + -(0-^r) 
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£=0 m=-l 
1 



P-1/2 
£m 



tm \ 



(34) 



1/2/ 



VT+b-p- 1/2' 
(_l)^2P+i/2 ^ 



n=0 



T{p+l/2)T{p + n + l/2) 



for p + i < 0, 
for p + i > 



T{p + n-l+ l/2)T{p + n + l + 3/2) n! V 1 + & ' 
We note that, although what we need here is only the case of an integer p, the above 
formula is valid for any p (except for the case p = —1/2). 

After the decomposition, we can take the radial coincidence limit r — > tq (followed 
by the angular coincidence limit if desired). Here we briefly explain the reason why the 
terms proportional to R/^ and 0/^ give no contribution to the final result. The term 
R/^ corresponds to R times the case of p = 0, for which -D^"^''^ is finite in the limit 
a — > 1 (i.e., R ^ 0). Hence all the coefficients vanish in the radial coincidence limit. As 
for 0/^, it can be replaced by (0 — 0')/^ which is equivalent to di^^ in the coincidence 
limit. This corresponds to the case of p = 1 multiplied by m. Hence all the harmonic 
coefficients become odd functions of m, and their sum over m for each i vanishes in the 
angular coincidence limit. As a result, the non- vanishing contribution comes only from 
the terms R/^\ 0/e^ 1/^, <^V^^ and 07^^. 

Barack and Ori define the regularization counter terms as 

A^L + B^ + CJL + 0(L-2) 



hm F„t 



D, 



E 

1=0 



lim F, 



dir 

al 



AfyL Bf, 



CJL 



(36) 
(37) 



where F^f 



is the multipole /-mode of F^^^ , L = i + 1/2, and Aa, and Ca are 
independent of L. The Aa term is to subract the quadratic divergence, the Ba term 
the linear divergence, and the term the logarithmic divergence. The term is the 
remaining finite contribution of the direct force to be subtracted. We find Ca = F>a = 
in agreement with Barack and Ori. We also find the complete agreement of Aa and 
Ba terms with their results for a general geodesic orbit as given below. We call these 
feature of regularization counter terms the standard form. 

The A-term describes the quadratic divergent terms of the direct force. It comes 
from the terms R/^^ in Eqs. fl30p ~ ()33|) . The important fact is that it is odd in R. 
This leads to the harmonic coefficients proportional to sign(i?). We find 

2M 



At = sign(i?) 



Ur 



2 ' 



A^ = 



sign(_R) 



ro 1 + Cyri 



r, 



ro - 2M 1 + C^/rl 



2 ' 



(38) 

(39) 
(40) 
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These A-terms vanish when averaged over both hmits R ±0. 

The 5-term describes hnearly divergent terms, which are of the form, 0^"/ ^^n+i 
Eqs. fl3()|) ~ fl33|) . We find the 5-term in terms of the hypergeometric functions as 







2r3 V2'2' ' rg 

1 /I 1 £2 /3 3 £2 



^U'U'o;i;--r +^^U'o'2;--t) ) , (42) 



2r^ V V2 2' ' r^y 2r^ V2 2' ' 
(ro-2M)£u.. / /3 3 
V l2'2'^'"7f 

3 3. £2X g£2 ^ /5 5 £ 



As mentioned at the beginning of this section, the above results for the A and B- 
terms perfectly agree with the results obtained by Barack and Ori in a quite different 
fashion [201 1^1122] • 

3. Gravitational Self- force in the Regge- Wheeler gauge 

Recently, Detweiler and Whiting found a slight but important modification of the direct 
and tail parts of the metric perturbation ^T]. The full metric perturbation in the 
harmonic gauge is now decomposed as 

^?r=^^J^+^^^^ (44) 

where the superscript H stands for the harmonic gauge, and the S part, h^J^, satisfies 
the same linearized Einstein equations with source as h^^l"^^-, 

□ fvl;^ + 2R,\Phlf = -WnT,, , (45) 

where hap = h^p — \ga(3h^^. The new tail part, called the R part, h^^, is then a 
homogeneous solution. Detweiler and Whiting showed that the properly regularized 
self-force is given by the R part of the metric perturbation. Namely, we have 

where z°'{t) is an orbit of the particle parametrized by the background proper time (i.e., 

g^u{dz^ / dT){dz^ / dr) = —1), and 

Fa[h] = -fiPj{hi3r,s - ^gp-yh^-s - ^h^s;i3 + ^g^sh\.p)u^u^ , (47) 

where Pa'^ = 5cf + UqM^ and = dz"" jdr. 

Since the R part is a solution of the Einstein equations, it is perfectly legitimate 
to consider a gauge transformation of it to another gauge. Therefore we can define the 
R part in an arbitrary gauge by a gauge transformation of the R part in the harmonic 
gauge to that gauge. We thus have an unambiguous definition of the self-force in an 
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arbitrary gauge. In particular, we may consider the self-force in the Regge- Wheeler 
(RW) gauge, in which we may be able to obtain the full metric perturbation. 

In this section, we formulate a method to obtain the self-force in the RW gauge 
by applying the mode-decomposition regularization. Then, as a simple but non-trivial 
example, we consider a circular orbit and derive the regularized self- force in the RW 
gauge to IPN order j^Hl- 



3.1. Gauge transformation to the RW gauge 

The gravitational self-force acting on the particle is given by the R part in the harmonic 
gauge, formally obtained as 

F^{t)= hni F^/^^'"] (x) 

= lim - h^-^] ix) 

= lim (F, {x) - [/i^'H] (x)) , (48) 

The S part can be calculated by the method discussed in the previous section. It is 
important to note that the S and R parts independently satisfy the harmonic gauge 
condition. 

We consider the gauge transformation of the R part from the harmonic gauge to 
the RW gauge, 

< ^ = + e^^""^ . (49) 
/^^"-^Y^=/^^"-2vef-r, (50) 

where is the generator of the gauge transformation. Then the self-force in the 

RW gauge is given by 

F^'^ir) = lim F„ [/i^'^wi 



lim F„[/i^'«-2Ve"^^W[/i^'«]](x) 

x— >z(r) 

lim F„ - /iS'H _ 2 V^H-^I^W n full,H _ ;^S,H1 1 

x— >2;(t) 

lim Fj/i™-2Ve"^^W[/i™l 
lim fF„[/if"™] 



(X) 



F„[/^S'H_2Ve^^^^[/i''«]](x)), (51) 



where we have omitted the spacetime indices of h^y and V (^^u) for notational simplicity. 
The full metric perturbation can be calculated by using the Regge- Wheeler- 

Zerilli formalism, while the S part h^^^ can be obtained with sufficient accuracy by the 
local analysis near the particle location. Thus provided that the gauge transformation 
of the S part can be unambiguously defined, we will be able to obtain the regularized 
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self-force in the RW gauge. Namely, the following gauge transformation should be well 
defined: 

Note that the self-force (|3T|) is almost identical to the expression obtained in the 
intermediate gauge approach if we replace the S and R parts by the direct and tail 
parts, respectively. The only difference is that the S and R parts are now solutions of the 
inhomogeneous and homogeneous Einstein equations, respectively. Hence the S part in 
the RW gauge is (at least formally) well-defined provided that the gauge transformation 
of the S part, Eq. (jH21) is unique. As will be shown later in Eqs. (IHUj). this turns out to 
be indeed the case. Therefore one may identify the self-force ()5H1 to be actually the one 
evaluated in the RW gauge [2H], not in some intermediate gauge. 



3.2. Full force for circular orbits 

We consider the full metric perturbation and its self- force in the case of a circular orbit. 
Hereafter, we will need to consider the case of general, non-circular orbits. However, 
because of several technical problems, which are not particular to the gravitational self- 
force but common to all types of the self-force, we leave it for future work. A promising 
method to deal with general, non-circular orbits will be presented later in Sec.|31 in which 
the case of scalar self-force is discussed for simplicity. An extension to the gravitational 
case is under progress |^ . 

First, the metric perturbation is calculated by the Regge- Wheeler- Zerilli formalism 
in which the metric is expanded in terms of Fourier series in frequency uj and tensor 
spherical harmonics with eigenvalue indices (£, m) by using the symmetry of the 
background spacetime. Then, we derive the Fourier-harmonic components of the self- 
force. As noted earlier, although the full self-force diverges on the orbit of the particle, 
each Fourier-harmonic component is finite. We may even sum over uj and m and the 
result is still finite. Thus we formally express the full self-force in the coincidence limit 
X z{t) as 

oo 

where r = tq is the radius of a circular orbit. 

On the Schwarzschild background, the metric perturbation h^^ can be expanded in 
terms of tensor harmonics as 

h = ^ f{r)Hairn{t, r)&f^ - iV2Hum{t, r)aJ2 + -^r^H2tm{t, r)sitm 
- ^vWTl)/^il(t,r)bH + iv/2£(m)/.Sl(t,r)b,„ 
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+ (^Ke.rn{t,r) - ^^^^^ Gimit.r)^ ge^ 

hoim{t, r)c}^ H hum{t,r)Cim 



(54) 



where a^^, a^^; ■ ■ ■ are the tensor harmonics introduced by Zerilh fT^. The energy- 
momentum tensor of a point particle takes the form, 

dr dr 

= ^]_u^u''^±^^^s^'\n - no{t)) . (55) 

The RW gauge is defined by the conditions on the metric perturbation as 

^RW ^ ;^(e)RW ^ ^(e)RW ^ ^'^^^ ^ q _ (gg) 

The Regge- Wheeler and Zerilli equations are obtained by plugging the metric 
perturbation ()54|) into the linearized Einstein equations and expand it in the Fourier 
series. The Regge- Wheeler- Zerilli formalism has been improved recently by Jhingan 
and Tanaka jHl]. Here, however, we use the original formalism developed by Regge and 
Wheeler, and by Zerilli. The tensor harmonics can be further classified into even and odd 
parities. Even parity modes are defined by the parity (—1)^ under the transformation 
{9, (p) ^ {n — 9 , (f) + n) , while odd parity modes are by the parity (—1)^"''"'^. 

For odd parity modes, we introduce a new radial function R£l^^\r) in terms of 
which the two radial functions of the metric perturbation are expressed as (also see jSH! 
for hfZJ 

>RW D(odd) .p-^x 

'hemui ~ _ 2]\^^ ^emu) ■ W ) 

The new radial function R^^^\r) satisfies the Regge- Wheeler equation, 

d"^ ^(odd) , r 2 T./ MD(odd)_ 87U r -2M 



d 2M 



X I -r^ 



^^[(1-— )Dw] + (r-2M)[(£-l)(£ + 2)]i/2g,„^J , (58) 
where r* = r + 2Mlog(r/2M — 1), and the potential Ve is given by 

. (l _ (!(i±i) _ . (59) 



The source term Qimu} vanishes in the case of a circular orbit and 

(£+l)(£-l)(£ + 2)" 



-1/2 ^^4>Y 



Dimuj{r) = — - — /i— ^ 

2 J 

x6{r-ro)mdeY;^{9o,(j)o), (60) 
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(61) 



ro - 3M ' ro V ro - 3M 

where = u^/ = a/ M/r^ which is the orbital frequency. The orbit is assumed to be 
on the equatorial plane without loss of generality. 

For even parity modes with the parity (—1)^, we introduce a new radial function 

(7) 

R\jn^{r) in terms of which the four radial functions of the metric perturbation are 
expressed as (also see 



f„^ fTRW o-RW „„J o-RW N 



/^RW 

Imuj 



A(A + ly + 3XMr + QM"' 



R 



(Z) 



+ 



r-2M d 



,(z) 



r2(Ar + 3M) ^ 
r{r-2M)~ i{r-2Mf 



-C 



Ar + 3M ^^^'"^ ' r(Ar + 3M)' 
where A = (£ — 1)(£ + 2)/2 and the local source terms are given by 

47rv/2 Mr (1) 
Ar + 3M uj ^"^^ 



+ 



r 



(62) 



Curnw 



^^-41 + \b,^^ - lQr^r[\i{i + m- 1)(^ + 2)]-V2f, 



imuj 



+ 



iuo r -2M 
167rir^ 



ir 



r-2M 



-B 



£mu) 



li{i + l){i-l){i + 2)]-'/'Femu.. 

(1) 



(63) 



r - 2M'2 

Here the harmonic coefficients of the source terms Aimuj, ^Imui ^"^^ Bimuj vanish in the 
circular case and 



B 



(0) 
£mu> 



11 



-1/2 



2M\ 1 



-6{r -ro)mY;^{eo, 0o) , 



+ !){£- + 2)' 



-1/2 



X 5{r - ro) (£(£ + 1) - 2m') F4(^o, <Po) • 
The new radial function R^^^{r) obeys the Zerilli equation, 



where 



and 



(Z) 



2M\ 2A2(A + ly + 6A2Mr2 + 18AMV + ISM^ 
r3(Ar + 3M)2 



(64) 
(65) 
(66) 



(z) 

imuj 



r -2M d 



ir 



— i- 



r dr 
(r - 2M)- 



(r - 2M)2 
_r(Ar + 3M) - 2M 

A(A + l)r2 + 3AMr + 6M2 ~ 



r(Ar + 3M)^ 



C2£muj 
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Ar2 - 3AMr - SM^ ~ 



(r - 2M) "^^^"^H • ^^^^ 

The Zerilli equation can be transformed to the Regge- Wheeler equation by the 
Chandrasekhar transformation pSj if desired. 

The homogeneous solutions of the Regge- Wheeler equation are discussed in detail 
by Mano et al. jHE] • By constructing the retarded Green function from the homogeneous 
solutions with appropriate boundary conditions, namely, the two independent solutions 
with the in-going and up-going wave boundary conditions, we can solve the Regge- 
Wheeler and Zerilli equations to obtain the full metric perturbation in the RW gauge. 
Here, we consider the radial functions up to the first post-Newtonian (IPN) order. 

The radial function for the odd part of the metric perturbation is obtained for 
r < ro as 



(2£ + l)^(£+l)(£ + 2) V^o. 
where Vt = u'^ /u^ (see j^H] for r > ro). For the even part, the radial function is obtained 
for r < ro (also see ^H] for r > ro) as 

^(Z) SfimTTuV / r (£2 -2£-l) Mr 



(2£+l)(£ + 2)(£+l)a;V ro (^-l)rg 
r^ ^ {f' - I + A)rQr \ ^ 



^ ' (2£+3)ro ^£(2£-l)(£-l) ' 



- YLid,. 0o) . (69) 



2(£4 + £3_6£2_4£_4)M 

The metric perturbation in the RW gauge is obtained from Eqs. ()57|1 and (j62|) . 

Next, we calculate the self-force. We apply the 'mode-decomposition regularization' 
method, in which the force is decomposed into harmonic modes and subtract the 
harmonic-decomposed S part mode by mode before the coincidence limit x —>■ z{t) 
is taken. Since the orbit under consideration is circular, the source term contains the 
factor 6{uj — niQ), and the frequency integral can be trivially performed. Hence we can 
calculate the harmonic coefficients of the full metric perturbation in the time-domain. 
This is a great advantage of the circular orbit case, since the S part can be given only 
in the time-domain. We also note that the ^-component of the force vanishes because 
of the symmetry, and F"^ = [(ro — 2M)/(rQf2)]F* for a circular orbit. 

Performing the summation over m, we find in the end. 



-^fuU.RWl^ — -^full,RW|£ — -^fuU.RW 



0, 



-^fulLRW 



-^fulLRW 



(£ + l)/i2 1 /i2(i2£3 + 25£2 + 4£-21)M 

7l ^2 4 (2£ + 3) (2£-l) 
I fi'^{12e + llf-10i+12)M 



rg 2 (2£-l)(2£ + 3)r3 " ^^^^ 

We see that the only non-vanishing component is the radial component as expected 
because there is no radiation reaction effect at IPN order. In the above, the index 
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(+) denotes that the coincidence hmit is taken from outside (r > tq) of the orbit, and 
(— ) from inside (r < tq) of the orbit, and the vertical bar with £, ■ ■ ■ , denotes the 
coefficient of the C. mode in the coincidence hmit. We note that the above result is valid 
for i >2. There are some complications with £ = and 1 modes, and they need to be 
treated separately. The details are discussed in [^HIEZI- Here we focus on the modes 
^>2. 

3.3. The S part in the harmonic gauge 

In this subsection, we present the S part of the metric perturbation and its self-force 
by using the local coordinate expansion. It may be noted that if we simply apply the 
scalar harmonic expansion to each component of the self-force, we obtain the exactly 
same regularization counter terms as in the case of a scalar charge discussed in Sec. |21 
However, here we consider the tensor harmonic expansion. So, each i component may 
be different from the scalar case. 

The S part of the metric perturbation in the harmonic gauge is given covariantly 

as 

[ a.^{x,z,ct)u'y{T,ct) 

+ 0{y'), (71) 

where hf,„ = the trace-reversed metric perturbation, z^et = -^(Tret); Tret 

and Tadv are the retarded proper time and advanced proper time, respectively, defined 
by the cross section of the past and future light cones of the field point x with the orbit. 
As in Sec. |21 g^a is the parallel displacement bi-vector, and y is the difference of the 
coordinates between x and Zq, y'^ = x^ — Zq. Details of the local expansion are given in 
Ref. JH]- The difference between the S part and the direct part appears in the terms 
of 0{y), i.e., the second term on the right-hand side of Eq. (|7T|) . It turns out, however, 
that this term does not contribute to the self-force. 

The local coordinate expansion of the S part is then performed in exactly the same 
manner as in Sec. |21 with a suitable extension of local quantities over the whole sphere. 
It is in principle possible to calculate the harmonic coefficients of the extended S part 
exactly. However, it is neither necessary nor meaningful because the extended S part 
is only approximate. In fact, corresponding to the fact that all the terms in positive 
powers of y vanish in the coincidence limit, it is known that all the terms of 0(1/L^) 
or higher, where L = i + 1/2, vanish when summed over £ ^H] in the harmonic gauge. 
It should be noted, however, this result is obtained by expanding the force in the scalar 
spherical harmonics. In our present analysis, we employ the tensor spherical harmonic 
expansion. So, the meaning of the index i is different. Nevertheless, the same is found 
to be true. Namely, by expanding the S part of the metric perturbation in the tensor 
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spherical harmonics, the S force in the harmonic gauge is found to have the form, 



i^'^L + B'^ + D'^, (72) 



where and are independent of £, and the ± denotes that the hmit to Tq is taken 
from the greater or smaller value of r, and 

jjtl ^ " J , J , (^o\ 

^ (L2 - 1)(L2 - 4) (L2 - 1)(L2 - 4)(L2 - 9) ^ ' ^ ^ 

Then, the summation of D'^ over £ (from £ = to oo) vanishes. For convenience, let us 
call this the standard form. As we shall see later, the standard form of the S force is 
found to persist also in the RW gauge, at least in the present case of a circular orbit to 
IPN order. 

For the moment, let us assume the standard form of the S force both in the 
harmonic gauge and the RW gauge. Then, we may focus our discussion on the divergent 
terms. When we calculate the S force in the RW gauge, we first transform the metric 
perturbation from the harmonic gauge to the RW gauge, and then take appropriate 
linear combinations of their first derivatives. We then find that the harmonic coefficients 
^2^m' ^olm^ ^Um^ differentiated two times, and is differentiated three 
times, while the rest are differentiated once, to obtain the S force. So, it is necessary 
and sufficient to perform the Taylor expansion of the harmonic coefficients up to 0{X^) 
for hl'^, /if^l^'" and hf^^^, and up to 0{X^) for G'f;^, and the rest up to 0{X), where 
X — T{= t — to) or R{= r — ro). To the accuracy mentioned above, the harmonic 
coefficients of the S part are found in the form, for example, 

ht'^it, r) = -nii [ £(2)(L2-l) + " " " J ^im%. <Po) , (74) 

where we have defined 

£(2) ^ + 1) ^ |^i,2 _ 1 ^ _ (75) 

The above example is the coefficient obtained by approaching the orbit from inside 
(r < ro). The other coefficients, as well as the coefficients obtained in the case 
of approaching to the orbit from outside (r > ro), can be found on the web page: 
http: / / www2.yukawa.kyoto-u.ac.jp/~misao /BHPC / . 

Now we consider the S force in the harmonic gauge. It is noted that the t, 9 and 
0-components of the S force vanish after summing over m modes. Using the formulas 
for summation over m 



m 



m 

the r-component of the S force is derived as 

^s,h|^ = ^s,h|^ = ^stn = 0> 
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-^S,H 



1^2(2ro-3M) l/x2(4ro-7M) 

— -p — 5 — — 

2 Tq 8 Tq 

/i^ M (172 - 14784 + 299) 
^ 128r3(L2 - 1)(L2 -4)(L2 -9) 

2 Tq 8 Tq L"' 

This is indeed of the standard form. It should be noted that the factor ^C^^) -which is 
present in the denominators before summing over m is cancelled by the same factor that 
arises from summation over m. If it were present in the final result, we would not be 
able to conclude that the summation of D'^ over I vanishes. We note that, apart from 
the fact that the denominator of the term takes the standard form, the numerical 
coefficients appearing in the numerator should not be taken rigorously. This is because 
our calculation is accurate only to 0{y^) of the S force, while the numerical coefficients 
depend on the 0{y) behavior of it. It is also noted that the 0(l/L)-terms are absent in 
the S force, implying the absence of logarithmic divergence. 

3-4- The S part in the Regge-Wheeler gauge 

Now, we transform the S part of the metric perturbation from the harmonic gauge to 
the RW gauge. The gauge transformation functions are given in the tensor-harmonic 
expansion form as 

^(odd) ^ Y.^^im^^'^it^ r) Uo,^d^Yern{e,<j)),smedeYern{e,A , 
^(even) ^ ^|M5«(t, r)r,„(^, 0) , Mf£-^^(t, r)r,^(^,0), 

im 

^2i™(^, r)deYUO, 0), M2^,'™(t, r)d,YUO, 0)} • (78) 

There are one degree of gauge freedom for the odd part and three for the even part. To 
satisfy the RW gauge condition (|^. we obtain the equations for the gauge functions 
that are found to be rather simple, and then we find 



A^r'^C.--) =5'S(*.''). (79) 



^, I . (80) 

We note that it is not necessary to calculate any integration with respect to t or r. It is 
also noted that the gauge functions are determined uniquely. This is because the RW 
gauge is a gauge in which there is no residual gauge freedom (for £ > 2). 
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Then, the S part of the metric perturbation in the RW gauge is expressed in terms 
of those in the harmonic gauge as follows. The odd parity components are found as 



S,RW/ 



h 



oim (t^ r) = hlfjt, r) + dtA 



S.H^RW 



it, r) , 



A' 



S.H^RW 

£m 



it, r) 



and the even parity components are found as 



+ 



S,RW 

2r 



it, r) = H^l^{t, r) 



r -2M 

S,RW 



S,H^RW 
Oim 

S,H, 



it, r) 



M{r -2M) 



M- 



S.H^RW 
Um 



it, r) 



Hum it, r) = H{^^{t, r 



+ 



dfM- 



S,H^RW 
Um 



(t, r) + dr M, 



S,H^RW 
Oim 



it, r) - 



2M 



r(r - 2M) 



m; 



S.H^RW 



h: 



+ 



S,RW 

2(r - 



it, r) 
2M) 



drM\ 



S,H^RW 
Um 



it, r) + ■ . 

r[r - 

2{r -2M) 



M 



2M)^^^™ 



S.H^RW 



Oim 



it, r) 



it, r) 



m: 



S,H^RW 
Um 



it, r) , 



(82) 



where the gauge functions are given by Eqs. (fTUj) and (jHUj) . 

Inserting the S part of the metric perturbation obtained as described in the previous 
subsection to Eqs. (f7^ and (jSUj) . we obtain the gauge functions that transform the S part 
from the harmonic gauge to the RW gauge. It may be noted that the gauge functions 
do not contribute to the metric at the Newtonian order. In other words, both the RW 
gauge and the harmonic gauge reduce to the same (Newtonian) gauge in the Newtonian 
limit. The S part of the metric perturbation in the RW gauge is now found in the form, 
for example, 



Km it,r) = j7rfi 



4:iTmro{L^ -2) {u'f'f 



de ^imi 



Next we calculate the S part of the self-force. Of course, it diverges in the 
coincidence limit. However, as we noted several times, in the mode-decomposition 
regularization in which the regularization is done for each harmonic mode (£-mode), 
the harmonic coefficients of the S part are finite. The calculation is straightforward. 
We find that the t, 9 and 0-components of the S force vanish after summing over m. 
Summing the above over m, the r-component of the S force is derived as 



p* I 

-^S,RW| 
-^S,RW 



S,RW 
S,H 



S,RW 







M) 



We now see that the S force in the RW gauge also has the standard form as in the case 
of the harmonic gauge and there is no 0(1/L) term. This is because the difference of 
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the S force between the harmonic and RW gauge arises from 2PN order as 

M(ro-2M)2 4!^'^(to,ro) 



(ro-3M) 
\2 ^ n^,^i 



3 M(ro-2M)^ g^C.^^^fa, rp) 
2 rg(ro-3Af) 

5.5. Regularized gravitational self-force to IPN order 



YUK 0o) ■ (85) 



In the previous two subsections, we calculated the full and S parts of the self-force in 
the RW gauge. Now we are ready to evaluate the regularized self-force. But there is 
one more issue to be discussed, namely, the treatment of the £ = and 1 modes. 

The full metric perturbation and its self-force are derived by the Regge-Wheeler- 
Zerilli formalism. This means they contain only the harmonic modes with £ > 2. If 
we could know the exact S part, then the knowledge of the modes £ > 2 would be 
sufficient to derive the regular, R part of the self-force, because the R part of the metric 
perturbation is known to satisfy the homogeneous Einstein equations [TT], and because 
there are no non-trivial homogeneous solutions in the £ = and 1 modes. To be more 
precise, apart of the gauge modes that are always present, the £ = homogeneous 
solution corresponds to a shift of the black hole mass and the I = \ odd parity to 
adding a small angular momentum to the black hole, both of which should be put to 
zero in the absence of an orbiting particle. As for the I = \ even mode, it is a pure 
gauge that corresponds to a dipolar shift of the coordinates. In other words, apart from 
possible gauge mode contributions, the £ = and 1 modes of the full force should be 
exactly cancelled by those of the S part. In reality, however, what we have in hand is 
only an approximate S part. In particular, its individual harmonic coefficients do not 
have physical meaning. Let us denote the harmonic coefficients of the approximate S 
force by F^''^^, while the exact S force and the full force by and respectively. 
Then, the R force F^ may be expressed as 

f^=j: {Fr - f!) = J2 (Fr - f!) 

i>2 e>o 
= j:{Ft^'-F!'^-)-J2D, 

i>2 £=0,1 

where Di = Ff — F^''^'^, and the last line follows from the fact that F^'^"^ are assumed 
to be obtained from a sufficiently accurate spherical extension of the local behavior of 
the S part to guarantee J2e>o F>£ = 0. Thus, it is necessary to evaluate the £ = and 1 
modes of the full force to evaluate the self-force correctly. 

First, we consider the contributions of £ > 2 to the self- force. As noted before, 
for the IPN calculation, the only r-component of the full and S part of the self-force 
is non-zero. The £ mode coefficients corresponding to the first term in the last line of 
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Eq. (|HH|) are derived as 

rpr I 
-^full,RW| 

45 /i^M 



I rpr rpr 

-^RWU ~ -^full,RW|^ -^S,RW| 



8(2£- l)(2£ + 3) rl ' 
Summing over £ modes, we obtain 

3/i2M 



(87) 



i^Rw(^ > 2) 



4rg 



Next, we consider the £ = and 1 modes. It is noted that the £ = and 1=1 
odd modes, which describe the perturbation in the total mass and angular momentum, 
respectively, of the system due to the presence of the particle, are determinable in the 
harmonic gauge, with the retarded boundary condition. On the other hand, we were 
unable to solve for the £ = 1 even mode in the harmonic gauge. Since it is locally a 
gauge mode describing a shift of the center of mass coordinates, this gives rise to an 
ambiguity in the final result of the self-force. Nevertheless, we were able to resolve this 
ambiguity at Newtonian order, and hence to obtain an unambiguous interpretation of 
the resulting self-force. 

The correction to the regularized self- force that arises from the £ = and 1 = 1 
modes, corresponding to the second term in the last line of Eq. ()86|) . is found as [231211 

,fJ„(,.0.1)^^-5Zi^. (89) 
'o ^'o 

Finally, adding Eqs. (jHH|) and (jEHl), we obtain the regularized gravitational self-force 
to the IPN order as 

'O 'O 

Since there will be no effect of the gravitational radiation at the IPN order, i.e., the t- 
and (/)-components are zero, the above force describes the correction to the radius of the 
orbit that deviates from the geodesic on the unperturbed background. It is noted that 
the first term proportional to /i^ is just the correction to the total mass of the system 
at the Newtonian order, where tq is interpreted as the distance from the center of mass 
of the system to the particle. Although this is a bit cumbersome way to obtain the 
Newtonian force, it is nevertheless an important result because it was derived from the 
relativistic self-force for the first time. 



4. New Regularization 

As mentioned in the previous section, the frequency spectrum becomes monochromatic 
in the circular case so that the frequency integral can be trivially performed to give the 
harmonic coefficients of the full metric perturbation in the time-domain. But when we 
consider the general orbit, the frequency spectrum is highly non-trivial and the frequency 
integral becomes formidable to perform analytically, making the S part subtraction 
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difficult to carry out. In this section, we present a new analytical method that can 
circumvent this difficulty |38j . 

Here, again for simplicity, we go back to the case of a particle with a scalar charge 
g, and consider the self-force given by Fjytp] defined in Eq. (jl)). We ffist describe a 
new decomposition, which we call the S-R decomposition, of the Green function in the 
Fourier-harmonic domain. Next we show that only the S part needs to be regularized, 
but it can be easily transformed to an expression in the time domain if we employ the 
post Newtonian expansion. We then calculate the (S— S) part, which is finite and regular, 
for general orbit to 6PN order and for circular orbit to 18PN order up to now, although 
we present the explicit result only to 4PN order in this paper. Finally, specializing to 
the case of circular orbits, for which we can calculate the R part analytically as well, we 
evaluate the self-force for several radii. In particular, we compare our result with the 
result numerically obtained by Detweiler, Messaritaki, Whiting and Diaz- Rivera [!?n|l4()j. 

4.1. The S-R decomposition 

The full scalar field induced by a scalar-charged particle is given by Eq. in terms of 
the retarded Green function G^'^^^{x,x'). The retarded Green function is represented in 
terms of the Fourier-harmonic decomposition as 

G^-^^, = / ^ e-(*-*') ^?^i^Ur, r')Y,^{e, mUQ' . 0') • (91) 

Then, the Klein-Gordon equation Q reduces to an ordinary differential equation for 
the radial Green function as 



^ 2M\ dP ^2{r - M) d ^ [ ru? £(£+1) 



(ir^ ' dr \r — 2M 



^5{r - r') . (92) 



The radial part of the full Green function can be expressed in terms of homogeneous 
solutions of Eq. (j^ . which can be obtained using a systematic analytic method 
developed in Ref. [3^1 • We have 



X «(r)<p(r')^(r' - r) + <p(rX(r')^(r - r')) 

W^w(<^rn,<^up) 



r^l 1 



2M 



(93) 



Here, the in-going and up-going homogeneous solutions are denoted, respectively, by 
ip\^ and </9^p, and v is called the 'renormalized angular momentum' ^^14111^ . which is 
equal to ^ in the limit Muj — > 0. 

We express the homogeneous solutions ip^^ and v?up terms of the Coulomb wave 
functions < and y^-'^-i gH US iSI ; 

^^ = ^:+p,^-^-\ ^i^ = %^:+^-''-\ (94) 
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The properties and the relations of the coefficients (3^, % are studied in Ref. [H] and jl21- 
An essential point to be noted here is that when we consider the post Newtonian 
expansion, i.e., when (f'^ is expanded in terms of z := ur and e := 2Muj, := {2z)~'^(p'^ 
contains only terms that are integer powers of z and e. Furthermore, this PN expansion 
turns out to be a double Taylor series expansion in and e/z, i.e., there appear only 
positive powers of u^. 

We now divide the Green function into two parts, as 



where 



aln^r, r') 



(95) 



-1 



X 



^:{r)^~''-\ry{r' - r) + ^-'^-\r)^:{r')e{r - r') , (96) 
-1 



X 



(^:(r)^-^-i(r') + v-'^-\r)^:{r')) 



+ %^:{r)^:{r') + /3.^,^^-^(r)^-'^-^(r') 



(97) 



Using results derived in Ref. jH], we obtain the behavior of the coefficients {I3y, •y^} in 
the PN expansion as 

P^ = 0{v''^'), ^, = 0{v-'), (98) 

where we have set z = 0{v) and e = 0{v^) with v being the characteristic orbital 
velocity. The functions ip'^ and ^Pc'^'^ are, respectively, of 0{v^) and 0{v~^~^) (except 
at £ = 0). Therefore, the three terms in the R part of the Green function become, 
respectively, of 0{v^^), 0(f^^~^) and 0{v^^~^^) relative to the S part. 

As should be clear from Eqs. ()9(ij) and ()97j) . the R part is regular at r = r', hence 
the part that needs to be regularized is the S part. Furthermore, if we truncate the PN 
expansion at a finite order, the R part terminates at finite i. Therefore the expression 
for the R part of the self-force in new decomposition takes following form. 



(A 



R 

a 



(99) 



£=0 



4-2. Computation of the S part 

We now compute the force due to the S part for a general orbit. The difference between 
the S force and the S force should be finite, because the R force is finite. Thus, in 
general, when expanded in terms of the spherical harmonics, the S force must take the 
same form as the S force 



(100) 
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Below we confirm explicitly that both and for the S force coincide with those of 
Eqs. (pn|) and ()43|) . Therefore, the S force minus the S force, which is finite, is given by 



;iS-S 



£=0 



5: ft 

1=0 



(101) 



To obtain an expression for the S force in the form of Eq. ()100|1 . it is necessary 
to perform the uj integration explicitly. Here, the key fact is that there appears no 
fractional power of uo in the S part. This is because we have chosen and 0^"^"^ as 
the two independent basis functions. As noted above, except for the overall fractional 
powers and z~^~^, they contain only the terms with positive integer powers of u"^. 
When we consider a product of these two functions, uj contained in the overall factors z" 
and z~"~^ just produces o;"^, which is canceled by uj from the inverse of the Wronskian. 
Thus, gemuj{r,r') is expanded as 



k=0 

The Fourier transform of cj^" simply produces 



(102) 



Differentiation of the delta function in the expression above can be integrated by parts 
to act on the source term. Thus, the integration over uj can be performed easily, and 
we can express the S force in the time domain as 



m, k 



g2 hm P„^V/3 

x— >z(t) 

{-if{d,r'^G,^,{r, z^itwuo, mj^'it), z^t)) 



. (103) 



The transformation of the S part into the time domain makes it possible to subtract 
the divergent S part analytically. If we were to perform this subtraction numerically, 
the fraction to be subtracted would become closer and closer to unity as i increases. 
Apparently, this would mean a stringent requirement on numerical accuracy. In this 
sense, we anticipate a clear advantage in the analytical subtraction. 



4.3. The (S-S) force 
The result for the S force is 



and 



2 ^. 00 

n=0 



S 



c 



q u^C 

U n=0 



(n) 



(104) 
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where the coefficients -f^^"), whose superscript {n) represents the PN order, are formally 
given by 



in) 



K 



in) 



E 1: 

i+3+k=n 
i+j+k=n 



(ijk) 



(ijk) 



(106) 



Here, the quantities cij*^*^'' are some functions of £, tq = z'^{to), 5s = 1 — U = — , 

C Tq 

and S and C are, respectively, the energy and the angular momentum of the particle. To 
obtain these expressions, we have used the first integrals of the geodesic equations and 
we have also reduced higher-order derivatives with respect to t by using the equations 
of motion. 

To summarize, the (S — S) force is given by 



?s-s 



E 

oo 

E 

i=0 



- F, 



F: 



-Br 



The final result is 



and 



nS-S 



nS-S 



" n=0 



Eci 



S-S(n) 



nS-S 



0, F 



s-s 



" n=0 



Ec^: 



S-S(n) 



_i?s-s _ 
ii'^(l-2M/ro) * 



C 



nS-S 



(107) 



(108) 



(109) 



where the coefficients Ca~^^"'\ whose upper index, n, represents the PN order. The 
calculation of Ca ^^"^ to a given order of post Newtonian expansion is straightforward. 
Most importantly, once we compute them, we can just keep the results and apply them 
to any general orbit at any time. Up to now, we have computed the (S — S) force to 6 
PN order. As example, we list some low PN orders terms: 

73 



^S-S(O) 



^S-S(l) 
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Once we obtain the general expression for the (S — S) part of the force, computation 
of the remaining R part is rather easy, because only terms up to a finite value of i 
contribute to the force for a given PN order. Then, the R force, which is what we want 
in the end, is given by the form (j^^ . 

One of the important properties of the (S — S) force is that it contains only the 
conservative part of the force. From Eqs. ()104|) . the £-mode of the S force takes the 
form 

F,s|^ = ^i,,(r,^,£)M^ f||^ = J-,,,(r,£,£), = ^^,,(r, ^, £) M^(110) 

The S part of the force is known to have exactly the same form. This implies that the 
(S — S) force also takes the same form. Thus, after summing over i, we conclude that 
the final form of the (S — S) force is 

F^s=^,(r,^,/:)w^ F^s = ^,(r,£,£), F^^ = J^4r,S,C)u\ (111) 

We can now explicitly show that the above form of the force implies the absence of 
a dissipative reaction effect. In other words, the force is conservative. The equations of 
motion to 0(/i^) are given by 
D 



dr 



(112) 



where is the perturbed four velocity and D/dr is the covariant derivative. Then, we 
obtain the evolution equation for the perturbed energy £ := —fit^u'^ as 



dS 

dr 



dr 



(113) 
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where i^^ = (dt)^ is the time-hke Kilhng vector. This equation is integrated to give 



S = S 



J-'t{r)dr. 



(114) 



Here, £ is an integration constant, which we can interpret as the unperturbed energy. 
In the same manner, for the perturbed angular momentum £, we obtain 



C = C + 



T^{r)dr. 



(115) 



Thus we find that there is no cumulative effect on the evolution of the energy and 
angular momentum of the particle. In other words, a force of the form (jlllj) preserves 
the presence of the constants of motion 8 and £. Concerning the radial motion, vJ' can 
be expressed in terms of fiut = S and ^u^ = Chy using the normalization condition of 
the four velocity. We have 



± 



^2_(l_2M/r) (l + C'^/r 



1/2 



(116) 



Thus, m'" is obtained as a function of r. 



4--4- Testing the efficiency for circular orbits 

To examine the efficiency of this new regularization method, we revisit the problem of 
the self-force for circular orbits. For this purpose, we first calculate the S force for a 
circular orbit to 18PN order. Then, combining with the calculation of the R part, which 
can be done with sufficient accuracy in order not to spoil the 18PN order accuracy of 
the (S — S) part, the regularized scalar self-force is evaluated and is compared with 
the result obtained by Detweiler, Messaritaki and Whiting [39 , and very recently by 
Diaz- Rivera et al. 

The components of the (S — S) self-force, F^~^ = — F^, have been given for 
general orbits in Eq. ()108p . Since the (S — S) force is expressed in terms of local quantities 
of the particle, i.e., its position and velocity, what we have to do is just to specify the 
orbit. In the present case, we consider a circular orbit, given by Eq. Then we 

find only the r-component is non-vanishing for the (S — S) force because the t- and 0- 
components are directly related to rates of change of the energy and angular momentum, 
and are purely dissipative for circular orbits. To 4PN order, it is given explicitly by 

^2 - 73 16151_, 395567, 
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where V = ^y M/rg = roil. 

The components of the R force are formally given by 
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,7r 



1 TT^R 



118) 



The R force in this case can be completely obtained analytically because the integration 
with respect to u is done by substituting mfl for uj. Also note that if we need the 
precision up to n-PN order inclusive, it is sufficient for us to calculate the modes up to 
i < n + 1. The 4PN results, after summation over £-modes, are 
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(119) 



Here 7 is Euler's constant, 7 = 0.57 ■ ■ ■. As expected, the t-component which represents 
the energy loss rate starts at 1.5PN order, corresponding to dipole radiation. 

An important property of the R force is that each £-mode of it may be evaluated 
without performing the post Newtonian expansion. In other words, we do not have to 
expand the R part of the Green function in powers of uo. Instead, we can easily compute 
each £-mode with sufficient accuracy for any given radius rg. Thus, as far as the R part 
is concerned, what we actually employ is not the standard post Newtonian expansion, 
but the truncation of the series expansion in £ at a given ^max compatible with the PN 
order of our interest. 

The total self- force is obtained by summing the (S — S) force and the R force. We 
have computed the force accurate to 18PN order so far. But here, we present the result 
explicitly only to 4PN order, 

4 7 o 4. 2~ 



i^R 



47rr^ 



+ 



-7 H vr 

3 ' 64 

/604 297r2 66 , , , 14, 14 

+ In 2 In 7 

V 45 1024 5 ^ ^ 3 3 ' 



(120) 
(121) 



In the r-component of the scalar self-force there is a significant cancellation between the 
(S — S) part and the R part, and the total force begins at 3PN order. 

The accuracy is actually limited by the (S — S) part. In Fig. ^ we show the 
convergence of the r-component of the (S — S) force as a function of the PN order for 
several representative orbital radii tq, so the accuracy of the full regularized force can be 
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Figure 1. The relative error of the post-Newtonian formulas in the r-component of 
the (S - S) force for the radii ro = 6M, lOM, 20M and 50M. The horizontal axis 
is the order of post-Newtonian expansion. The top figure shows the convergence in 
the Taylor expansion and the bottom figure is the one obtained by using the Pade 
approximation. 



read from this figure. Here an estimator of convergence of the PN expansion is defined 

by 



irS-si T?s-s\ 

\{n)--t'a l(n-l) 



(122) 



where -Fa~^|(n) denotes the (S — S) part of the force truncated at n-PN order (inclusive), 
and the denominator Fa denotes the exact (fully relativistic) self-force including the R 
part. In practice, since it is impossible to know the exact value of it, we use the most 
accurate result in our calculation. It is found that the convergence of PN expansion 
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is steady even near ISCO, although it slows down there. The convergence improves 
slightly by using the Pade approximation near ISCO. Here, in the Pade approximation, 
we have chosen the denominator to be quadratic in V^. 

Now let us compare our result with the one obtained by Detweiler and his 
collaborators |Hn|l4()j. They calculated the radial component of the self-force for various 
radii in units of M = 1 and = Att. In our case, we employed a Pade approximation 
for the (S — S) force accurate to 18PN order and used the most accurate R force in 
our calculation including the terms up to £ = 19. Their results are compared with ours 
for tq = 6M (ISCO), lOM and 20M in Tabled As clear from it, the agreement is 
impressive for ro = lOM and 20M, while there is relative error of ~ 10~^ for tq = 6M. 
This is consistent with the error estimate given in Fig. ^ This error may seem large, 
but if we use our result as a template for a space gravitational wave detector such as 
LISA, it turns out that the error is small enough (44^. Thus, we conclude that our new 
method is capable of computing the self-force with sufficient accuracy, and the result 
obtained to 18PN order seems accurate enough even in the limit of ISCO. 



'"0 


6M 


lOM 


20M 




1.676820878 x 10"^ 


1.378448171 x 10"^ 


4.937905866 x 10"^ 


Detweiler et al. 


1.6772834 x 10^^ 


1.37844828(2) x 10"^ 


4.937906 X 10"^ 



Table 1. The r-component of the self- force in units of AI — 1 and — An obtained 
by our method (top column) and by Detweiler et al. [391 QUI (bottom column) . The 
number for vq — lOM is taken from Ref. |39| in which they estimate the error by 
a Monte Carlo simulation, and the numbers for ro = 6M and 20M are taken from 
Table I of Ref. 



5. Conclusion 

In this paper, we have reviewed recent progress in the regularization and computation 
of the self-force acting on a particle orbiting a black hole. 

We have presented a method to regularize the self-force analytically, employing the 
post Newtonian expansion. If we were to perform regularization numerically, the fraction 
to be subtracted would become closer and closer to unity as i increases. Apparently, this 
would mean a stringent requirement on numerical accuracy. This is a clear advantage 
in the analytical approach over a numerical one. Thus, as far as the Schwarzschild 
background is concerned, we are almost ready to calculate the gravitational self-force for 
general orbits with a sufficiently high accuracy in post Newtonian expansion. Research 
along this direction is now in progress |45i . 

At the same time, however, we should mention that our analytic method will become 
ineffective in the very high frequency regime. Therefore, it seems important to develop 
an alternative numerical method that may play a complimentary role. 
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Furthermore, although these recent developments are substantial, they are not quite 
good enough because our final goal is to compute the regularized self-force in the Kerr 
background. As mentioned in the introduction, there is some progress in the case of 
the Kerr background [2S1 12Z] • However, we should admit that we are still at a primitive 
stage in the case of the Kerr background. 

Recently, extending the adiabatic orbital evolution in a more systematic manner, 
Mino proposed an alternative, possibly much more powerful method to deal with the 
gravitational reaction force to the orbit of a particle |42J. Perhaps, we should test his 
idea by applying it to orbits in the Schwarzschild background and comparing the orbital 
evolution with the one obtained from the self-force regularization method. Together 
with computation of the gravitational self-force for general orbits in the Schwarzschild 
background, we hope to come back to this issue in the near future. 
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